function fea_sol=inti(G,V,g,v,K)
     fea_sol=randi(K,1,length(g));
     G_cur=zeros(1,K);
     V_cur=zeros(1,K);
for j=1:length(g)
    while(1)
        fea_sol(1)=randi(K);
        car_no=fea_sol(j);
        if G_cur(car_no)+g(j)>G(car_no)  || V_cur(car_no)+v(j)>V(car_no)
            fea_sol(j)=randi(K);
        else
            G_cur(car_no)=G_cur(car_no)+g(j);
            V_cur(car_no)=V_cur(car_no)+v(j);
            break
        end
    end
end
end